Trionic phase of ultracold fermions in an optical lattice: A variational study 
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To investigate ultracold fermionic atoms of three internal states (colors) in an optical lattice, 
subject to strong attractive interaction, we study the attractive three-color Hubbard model in infinite 
dimensions by using a variational approach. We find a quantum phase transition between a weak- 
coupling superconducting phase and a strong-coupling trionic phase where groups of three atoms are 
bound to a composite fermion. We show how the Gutzwiller variational theory can be reformulated 
in terms of an effective field theory with three-body interactions and how this effective field theory 
can be solved exactly in infinite dimensions by using the methods of dynamical mean field theory. 
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I. INTRODUCTION 

The Bose-Einstcin condensation (BEG) was realized in 
atomic traps in 1995 by cooling down ^^Rb atoms in a 
magnetic trap^. Since then, experiments have been per- 
formed on a variety of ultracold alkali atoms ranging from 
Li to Csjii^i^i^i^i^ In addition to displaying new phenom- 
ena such as the BEC-BCS crossoveri^ and the bosonic 
Mott-transition^, these systems provide also clean and 
flexible realizations for basic theoretical models such as 
the Hubbard model i^ii^ 

Alkali metal isotopes with odd (even) number of neu- 
trons behave as fermions (bosons).— Although initial 
experiments were mainly done on bosonic systems, the 
degenerate Fermi systems have also been realized and 
confined to optical lattices in recent years due to ad- 
vanced sympathetic cooling techniques J^i^^iii^ Because of 
the Pauli exclusion principle, these fermionic atoms can 
display a variety of interesting phenomena and phases 
that do not have an analog in bosonic systemsii^ In the 
remainder of the paper we shall focus on fermionic sys- 
tems of three internal degrees of freedom and show how a 
quantum phase transition appears in this system, which 
is a simplified version of the color superconductor-baryon 
phase transition in quantum chromodynamics (QCD). 

Typical hyperfine couplings are larger than the stan- 
dard experimental temperatures used to study ultracold 
gases. In the absence of an external magnetic field, the 
hyperfine coupling aligns nuclear (I) and electronic (S) 
spin antiparallel to each other, and the hyperfine spin 
F = S -I- I is conserved. In finite magnetic fields, only 
the hyperfine spin along the external field is a good 
quantum number. The hyperfine spin thus provides 
an internal quantum number that we shall refer to as 
"color" henceforth. 

Systems with three internal quantum numbers are of 
special interest, since they are rarely observed in solid 
state physics. Such a three-component fermionic system 
may be created, e.g., by trapping the lowest three hy- 
perfine levels of ^Li atoms in all-optical setups in large 



magnetic fields (see Fig. [T]). 

Such three-component systems with weak interactions 
have been studied first in Ref. [isl . where it has been 
shown that for small attractive interactions, a color su- 
perfluid state emerges. This work has been generalized 
to incorporate three-body correlations at large attractive 
interaction strengths in Ref. [l6l. At the same time, results 
for the three-fermion problem in a single parabolic well 
and a mean field calculation to describe a two-component 
BEC-BCS crossover through a Feshbach resonance in ^Li 
appeared ] ^^1^^ 

In this paper, we shall study three-color fermionic sys- 
tems in an optical lattice. Optical lattices are realized by 
standing light waves, which create a periodic potential 
for the trapped atoms If the amplitude of the lasers is 
strong enough, then the atoms are localized to the min- 
ima of this potential and at low temperatures can only 
move by tunneling between the lowest lying states within 
each such minimum. 
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FIG. 1: (Color online) Sketch of the lowest lying hyperfine 
levels of ®Li in external static magnetic fields. The high-field 
seeker states have negative tangent, which cannot be trapped 
in magnetic traps. However, using all-optical setups, one can 
trap the hyperfine states denoted as |2), and |3), and thus 
creating a three-component quantum degenerate fermionic 
system. 



The trapped atoms, however, also interact with each 
other. The dominant interactions between alkah atoms 
are short-ranged electrostatic van der Waals forces, which 
can be well approximated by S potentials at low energy 
scales. If the scattering length of the interacting atoms 
is smaller than the lattice constant of the optical lattice, 
then these interactions arc restricted to a given lattice 
site. Correspondingly, the interacting system of atoms 
in an optical lattice can be accurately described by the 
following Hamiltonian: 



with the creation operator of a fcrmionic atom of 
color a = 1, 2, 3 at site i, and fiia — cf^Cia- In the tun- 
neling term, implies restriction to nearest neighbor 
sites, and the tunneling matrix element is approximately 
given by t = En{2l ^)s^l^e-^'^'\ where = ^ is 
the recoil energy, q is the wavevector of the lasers, m is 
the mass of the atoms, s = Vq/Eii, and Vq is the depth 
of the periodic potential.— We neglect the effects of the 
confining potential in Eq. ([1]), which would correspond to 
a site-dependent potential term in the Hamiltonian. The 
interaction strength Uap between colors a and (3 is re- 
lated to the corresponding s-wave scattering length, Oq^, 
as Uap = Efi Uapq y/sjirs'^^'^ Note that fermions 
with identical colors do not interact with each other. 

For the sake of simplicity, we shall first consider the 
attractive case with Ua/s = U < 0. This case could be 
realized by loading the ^Li atoms into an optical trap 
in a large magnetic field, where the scattering lengths 
become large and negative a^^ « « —2500 Qq, for all 
three scattering channels, 12, 13, and 23i^ 

Introducing the usual Gell-Mann matrices, X'^p {a = 
1,...,8), it is easy to see that global SU(3) transfor- 
mations exp{iJ2iJ2aai3 4'aC'l^X'^i3Cii3) commutc with the 
Hamiltonian, which thus also conserves the total number 
of fermions for each color, Na — J2i ^ia- This conserva- 
tion of particles is only approximate because in reality the 
number of the atoms in the trap continuously decreases 
due to different scattering processes. Here, however, we 
shall neglect this slow loss of atoms and keep the density 
Pa of atoms for color a as well as the overall filling factor 
P = I Eq Pa fixed. 

Let us first focus on the case of equal densities, pa = p. 
For small attractive U < 0, the ground state is a color 
superfluidfi^ atoms from two of the colors form the 
Cooper pairs and an s-wave superfiuid, while the third 
color remains an unpaired Fermi liquid. However, as 
we discussed in Ref. for large attractive interactions, 
this superfiuid state becomes unstable, and instead of 
Cooper pairs, it is more likely to form three- atom bound 
states, the so-called "trions". These trions are color sin- 
glet fermions, and for large \U\ they have a hopping am- 
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FIG. 2: (Color online) The ground states for \U\ < \Uc\ and 
\U\ S> \Uc\ can be calculated by perturbation theory. The 
former is a BCS-state which breaks the SU(3) symmetry, and 
the latter is a trionic state with three-particle singlet bound 
states. 

plitude. 



Furthermore, one can easily see that if two trions sit on 
neighboring lattice sites, then they increase the energy of 
each other by an amount V ^ t'^ /{2\U\). This is because 
the energy of an individual trion is decreased by quan- 
tum fiuctuations where one of the atoms virtually hops 
to one of the neighboring sites. These quantum fiuctua- 
tions are reduced if the two trions sit next to each other. 
Therefore, trions will tend to form a Fermi-liquid in any 
finite dimensions. This Fermi liquid state may be further 
decorated by charge density wave order at large values of 
\U\. Also, the Fermi liquid scale TpL of the trionic Fermi 
liquid should depend on the value of U, and at the tran- 
sition point, U = Uc, we expect it to go zero (see Figs. [5] 
and El). 

In order to get analytic expressions, we shall study the 
ground state in d = oo dimensions. Then, to reach a 
meaningful limit and to get finite kinetic energy, one has 
to scale the hopping as t = with t* fixed. In this 

limit, however, trions become immobile. Therefore the 
d ^ oo trionic states are well approximated as 

\Ta)^Y[c+c^c+\0), (3) 

where A denotes a subset of sites where trions sit. We can 
calculate the energy of this state in infinite dimensions: 
a single trion has an energy 3U, thus the energy of such 
a state per lattice site is given by Et/N = 3Up, with 
E-T the total energy of the system and N the number of 
lattice sites. 

Clearly, the two ground states obtained by the per- 
turbativc expansions have different symmetries: the su- 
perfiuid state breaks SU(3) invariancc, while the trionic 
state does not. Therefore, there must be a phase tran- 
sition between them. Note that, relying on symmetries 
only, this argument is very robust and carries over to any 
dimensions. In infinite dimensions, we find that trions 
arc immobile. However, this is only an artifact of infinite 
dimensions and in finite dimensions, a superfiuid-Fcrmi 
liquid phase transition should occur. 
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One could envision that some other order parameter 
also emerges and masks the phase transition discussed 
here. Preliminary results (not discussed here) suggest 
that indeed a charge density state forms at large values 
of \U\, but except for half filling, which is a special case 
not discussed here, we do not see any other relevant or- 
der parameter that could intrude as a new phase. One 
could, in principle, also imagine a phase with simultane- 
ous trionic and fermionic Fermi surfaces, similar to the 
one of Ref. [2l|. However, in contrast to the calculations 
of Ref. [U (which do not take the local constraints on the 
lattice into account), in our scenario, there is no Fermi 
surface at the quantum critical point. 

From the above argument, it is unclear whether the 
phase transition is of first or of second order. To address 
this question in infinite dimensions and give some more 
quantitative estimates for the relevant condensation en- 
ergies and the critical value of U /t, we need an approach 
that is able to describe the superfluid state and, at the 
same time, also accounts for three-body correlations. In 
the present paper, we construct a variational Gutzwiller 
wave function that is able to capture both correlations 
simultaneously. We then show how averages can be eval- 
uated by constructing an effective action that contains 
three-body correlations and how calculations can be ana- 
lytically done in infinite dimensions by using the method- 
ology of dynamical mean field theoryi^ Our method, 
which uses a single suitably chosen Gutzwiller correlator 
and exploits the cavity method of dynamical mean field 
theory, is in its present form less flexible than the multi- 
band method reviewed by Biincmann et. al.^ but it in- 
corporates three-body correlations in a very transparent 
way. Our analysis shows that, within the Gutzwiller ap- 
proach, the phase transition is of second order in infinite 
dimensions. We therefore expect it to remain of second 
order in any finite dimension above the lower critical di- 
mension. Nevertheless, as we discuss later, the quantum 
phase transition may become of slightly first order due 
to a not perfectly SU(3)-symmetrical interaction. 

Interestingly, our results also suggest that there is a 
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FIG. 3: (Color online) Schematic phase diagram of the at- 
tractive SU(3) Hubbard model. The color superfluid phase 
is stable below a critical temperature Tc, where the appear- 
ance of superfluid domains is expected. At large interaction 
strengths |?7| > \Uc\, color singlet trions emerge instead of 
the Cooper pairs and form a Fermi liquid at low tempera- 
tures T <Tfl. 



tendency to create an imbalance of the densities pa in 
the superfluid phase. The physical reason is simple: one 
can gain condensation energy by transferring atoms from 
the unpaired channel to the paired ones and thereby cre- 
ate ferromagnetic order as a secondary order parame- 
ter. For equal densities Pa = p and color conserva- 
tion, his can only happen if the atoms are segregated 
and domains are formed, as shown in Fig. [31 where we 
sketched the schematic phase diagram of the attractive 
SU(3) model away from half-filling. This picture, which 
is first proposed in Ref. [l6l . has also been confirmed by 
the Ginzburg-Landau-type effective field theoretical ap- 
proach of Ref. m. 

So far, we discussed the fully SU(3) symmetrical case. 
As we also demonstrate later, the phase transition dis- 
cussed above is not sensitive to having perfect SU(3) sym- 
metry. On the other hand, to form the superfluid phase, 
it is important to have approximately the same Fermi 
momenta for at least two colors. 

The rest of the paper is organized as follows. In Sec. [TTl 
we introduce the Gutzwiller ansatz for the ground state. 
In Sec. mil we reformulate the Gutzwiller expectation 
values as an effective path integral. In Sec. IIVI we derive 
a local action, which can be used to solve the effective ac- 
tion in the d ~ oo dimensional limit. In Sec. |Vl we first 
summarize the results for the SU(3) symmetric case and 
then generalize the approach by analyzing the Hamilto- 
nian with nonuniform interaction strengths in order to 
describe a system of ^Li atoms. In Sec. I VII we present 
a brief discussion of analogies with QCD. In Sec. IVIIl 
we present our conclusions. Some of the technical details 
can be found in the Appendices. 

II. GUTZWILLER ANSATZ 

To capture the color superfluid - trion transition, we 
approximate the ground state of the infinite-dimensional 
system by the following Gutzwiller-correlated wave func- 
tion: 

IG) = [][1 + (.9 - l)nan^2n^3] |BCS) . (4) 
i 

Here 5 is a Gutzwiller variational parameter, which in- 
creases (or decreases) the amplitude of terms in the 
"uncorrelated-state" which have triple occupancies. In 
the g ^ 00 limit, \G) becomes a superposition of trionic 
states similar to Eq. ([3]). We choose the uncorrelated 
ground state as a BCS-likc state with colors "1" and "2" 
forming Cooper pairs^^ 

|BCS)= n c+HK' +1^1^5+151^,2)10) ■ (5) 

k: ek<M3 k' 

The operators 5j^^ in Eq. (O diagonalize the first part of 
the Hamiltonian ([1]), and create fermions with momen- 
tum k, color a, and energy et = —tJ2a^^^' '^^^'^ 
vector a running over nearest neighbor sites. The coeffi- 
cients u1 = 5(1 -I- ^k/ \/Ck~t-~A2) and = i/l - ^'"'^ 
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the usual BCS coherence factors, with = Ck — A*i2- The 
"chemical potentials" /ia and /ii2 appearing in the wave 
function should be considered merely as parameters that 
are adjusted to fix the densities pa at a given value. 

To perform a variational calculation, we need to evalu- 
ate the Gutzwiller expectation value of Hamiltonian llj. 



{H)g = {G\H\G)/{G\G) , 



(6) 



at a given filling factor p, and then minimize it with re- 
spect to A, and eventually the density of the third 
color, Pa. 

We note that our variational wave function smoothly 
interpolates between the color superfluid {g ~ 1) and 
the trionic state {g oo). We can also compare the 
Gutzwiller energy to certain reference states, such as the 
Fermi sea {g = 1. A = 0), the Gutzwiller correlated 
Fermi sea 7^ 1, A = 0), or the uncorrelated BCS state 
{g = 1, A ^ 0), and thereby estimate correlation or 
condensation energies. 



III. EFFECTIVE GRASSMANN THEORY 

We first derive an effective action that can be used to 
replace the Gutzwiller expectation value of an operator 
O, 



{d)G = {G\d\G)/{G\G), 



(7) 



by a combination of the Grassmann path integrals. For 
this purpose, let us first rewrite the denominator of Eq. 
dZD as 

(G|G) = (BCS| + (52 _ l)n,in,2n.3]|BCS). (8) 



Here, the operator in the middle is in a normal ordered 
form (i.e., all cf^ appear to the left of the operators Cia)- 
Since the BCS state is a non-interacting state (in terms 
of bogoliubons), we can use Wick's theorem to evaluate 
expectation values appearing in Eq. Due to the nor- 
mal ordering, every cf^ and Cia appears in Eq. ([S]) only 
once. Therefore, wc can calculate the quantum mechan- 
ical expectation values as the Grassmann path integrals, 
where we replace normal ordered operators by the Grass- 
man variables as Cia i^ia and fjia and evaluate 
expectation values with an appropriately chosen action 

^0, 



)bcs {■ ■ ■)s„ = S'D-qVi^ 
^0 = -5 I] 



-5o 



(9) 



Here, wc introduced the "Nambu spinors" , '^i = 
iVia, Via)i ^ud thc propagator D'^j must be chosen so that 
it satisfies the conditions, 



BCS 
{CjfjCia) BCS 



^ ijafi ' 



with G° and F° the normal and anomalous Green's func- 
tions. It is easier to express in the Fourier space, 
where it is just a 6 x 6 dimensional matrix in the Nambu 
space that can be expressed as 



i^"(k) = 
We can thus rewrite (G|G) as 



GO(k) FO(k) \ 
F"+(k) -GO+(k) ) 



(G|G) = (n 



i+(5^-i)n 



(10) 



(11) 



So 



Note that in the above procedure, we have not doubled 
the Hilbert space since the integration is over the fields, 
f] and 77, and therefore the 5* and 'I' do not represent 
independent Grassman variables. 

Note also that the Green's functions G°^o = 



)bcs and 



ijafS 



{cjfiCia)BCS can easily be de- 
termined in terms of the BCS coherence factors. Fur- 
thermore, neither the Green's functions nor the Grass- 
man fields r]j0 have a time argument. These Grassman 
fields should thus not be confused with the Grassman 
fields that appear in the path integral expression of the 
density matrix or the time evolution operator. 

To make further progress, let us define the "triple oc- 
cupancy" as U = niini2n.a = Ha ViaVia- Using the prop- 
erties of the Grassmann variables, this quantity can be 
expressed as a function of the Grassmann vectors defined 
above 



4 (*^r3*^)^ 



(12) 



where T3 = a3'S>Saf3, and the third Pauli matrix 173 acts in 
Nambu space. Clearly, <| = 0. Therefore, exponentiating 
the product in Eq. (fTT|) . we obtain the relation 



(G|G) 



(13) 



where the full (interacting) effective action is given by 
S - -^Y.'^^[D"-%^J [g' l)Y.U . (14) 



Interestingly, the Gutzwiller variational parameter ap- 
pears in this effective action as a three-body interaction. 
For the uncorrelated wave function, (7=1, thc inter- 
action term vanishes. For trionic correlations one has 
g > 1, i.e., the effective three-body interaction is attrac- 
tive, while for g < 1 the effective interaction is repulsive. 

This procedure can be repeated with certain operator's 
quantum mechanical expectation values after proper nor- 
mal ordering. For local operators Oi, the Gutzwiller ex- 
pectation value is 



(G|d,|G)/(G|G) = 
(BCS|:d, :g 



(5'-l)U|BCS>/(G|G) , 
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where the GutzwiUcr normal ordered operator : Oi :g is 
defined as foUows: 



^^ :g=: [1 + (ff - l%]d,[l + {g- 1%] 



(15) 



and ii = niini2ha. Here, the numerator can also be 
converted into expectation values evaluated with S: We 
first use Wick's theorem and replace all operators by the 
Grassman variables and the average by (. . .)so- Then, we 
can insert the missing term [1 + [g^ — l)ti\ and reexpo- 
nentiate the product to have an average with the action 
S. 

With this procedure, the particle density and double 
occupation on a site can be written as follows: 

Ni = {fill + fia + na)G 

= (na + 7ii2 + "13)5 + 3(g^ - (16) 

Di = {hiini2 + ni2ni-i + niifii^)G 

= {■niini2 + ni2ni3 + niinja)^ 



+3(.g'-l)(t,)< 



(17) 



To evaluate the expectation value of the kinetic energy, 
we also need to evaluate expectation values of the type 
{cf^Cja)G for « 7^ j- We therefore define the density su- 
pcrmatrix (formulas here are valid for i ^ j), 



Pij Qi 
Qij —P.i 



+ ' 



(18) 



with the 3x3 matrices Pij and Qij defined as 



Pi^3,al3 = {ctj-3p)G = {ViaVjl3)s 

+ (1 - g)'^{ViaVjp{diadjp)))s (19) 



and 



Qi^j,af3 = {CiaCjf})G = {'>lia'>ljp)s 

-(1 - g){riiarijf3{dia + djfj))s 
+ (1 - g)'^{nia'njpdiadji3)s- 



(20) 



In these expressions, we defined the Grassman factors 
dia = Y[f3^a "-^/s ■ Clearly, for g ^ I, the matrix V 
is simply related to the unperturbed Green's function, 
'P{g = 1) = [-0°]^. As we shall now see, although the 
above expressions look rather complicated and contain 
formally three-body correlation functions in terms of the 
effective field theory, we shall still be able to express V 
in terms of the full single particle Green's function D, 
defined as 



A, = 



(1 F- 



(21) 



and its proper self-energy E. 

At this point, it is useful to define the generating func- 
tional. 



In 



PryPryexp ( -5 + X]^'^3*» 



(22) 



where the "current" Ii is a six-component Grassmann 
vector, i, = {Jii,Ji2,Ja,Jii,Ji2,Jiz), and make use of 
standard field theoretical methods.— As usual, the com- 
ponents of the Green's function can be obtained from F 
by functional derivation. 



G 



Pij a j3 



7-0 



(23) 



1=0 



The (improper) self energy S is defined by cutting off the 
bare lines from the dressed propagator: 



D, 



J-'ipJpqJ-^qj^ 



(24) 



while the proper (one particle irreducible) self-energy 
obeys the Dyson's equation, 



Sij — -t- ^ T^ipDpgSqj . 



(25) 



pq 



It is useful to visualize these self-energies in terms of the 
Feynman diagrams of the effective field theory, as shown 
in Fig. m although we shall sum up these diagrams later 
by using nondiagrammatic methods. As mentioned ear- 
lier, the Gutzwiller correlator generated a three-body in- 
teraction in the effective field theory, while anomalous 
Green's functions appear due to the superfiuid correla- 
tions. 

To generate the necessary identities to evaluate the 
expectation values in Eqs. (|19p and (|20p . we shift the 
Grassmann fields fji and 77, in the generating functional 
F as 



^^P'ipTsIp, 



(26) 



where A is a real valued parameter. Note that this equa- 
tion also determines how the field ^i is changed, since the 



GO. 



GO 





FIG. 4: Examples of first and second order self-energy di- 
agrams. Continuous lines denote the unperturbed normal 
Green's functions, the staggered lines the anomalous ones, 
while vertices stand for the three-body interaction, ~ {g^ — 1). 
In infinite dimensions, we shall sum up these diagrams by us- 
ing a nondiagrammatic method. 
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components of i and ^'^ are related. Proceeding with 
the functional derivation after this change of integration 
variables leads to an expression of the dressed propaga- 
tor D in terms of expectation values of the Grassmann 
fields related to the expectation values in Eqs. and 
([20|l (for details, see Appendix 1X1 ) . and finally compar- 
ison of these with the definition of the self-energy gives 
the following identities {i ^ j): 

- {g'-l)\3{^dfdf^,)sr3^S,,, (27) 

(g2 - l)(vl/,df *,)5r3 = ^l^P^^ (28) 

p 

p 

where df — dn + di2 + da = \{^iT3^iY. Using these 
identities, it is possible to express the density superma- 
trix V in terms of the improper self-energy. This relation 
is more transparent in the Fourier space, where it reads 

V\\^) = i5"(k) + C 

+ K)-^)««(^°«-jTt)- 

(30) 

where the superscript t refers now to the transposed ma- 
trix in color space, and C is a k-indepcndent 6x6 matrix 
related to the i = j contributions. This term does not 
play a role in the evaluation of the kinetic energy, 

^^ = E^k^^""('^)' (31) 
k a 

since we assumed only nearest neighbor hopping, and 
therefore J2k,a CaaCk = 0. 

IV. CAVITY FUNCTIONAL AND 
CONNECTION TO DYNAMICAL MEAN FIELD 
THEORY 

In the previous section, we showed that the Gutzwiller 
expectation values can be expressed by path integrals, 
which still need to be computed. These path integrals 
cannot be exactly evaluated in general. However, an im- 
portant simplification occurs in the limit of infinite di- 
mensions, d oo. There the off-diagonal {i ^ j) Green's 
hmctions decay as Dij ~ where \\i — j\\ de- 

notes the minimum number of steps that are needed to 
reach the lattice point i from the other lattice point 
Therefore, the off-diagonal {i ^ j) components of the 
proper self energy rapidly go to zero, ~ l/d'^ll*"^!!/'^, 
and vanish in the infinite-dimensional limit, where the 
self-energy becomes local, 

~ SijY!i . (32) 

Here, we already used the discrete translational symme- 
try of the lattice: S(i) = S. 



As E is a central quantity of interacting systems and 
allows one to compute virtually any other quantity, our 
main goal shall be to determine E. This step is done sim- 
ilar to the cavity method also used in dynamical mean 
field theory to extract local properties at a mean field 
leveli^^ In this approach one integrates out the Grass- 
mann variables on all sites except for the origin, 

i J PS^P'rye-^ = ^e-^^ , (33) 

where the prime indicates that one should not integrate 
over the variables tjq and 770 at the origin. Thus, by con- 
struction, the local action depends only on the variables 
(770,770), Sl = [770 , , and all local expectation values 
and correlation functions are invariant, 

{Oo)s = {Oo)s, ■ (34) 

Importantly, the local Green's function is also invariant. 

Since the hopping is scaled as ~ 1/ Vd, it can be shown 
(Appendix |B]) that in infinite dimensions, the generated 
local action takes on a simple functional form, 

-Sl^ l^PV^'^ + ^^i^T.^f , (35) 

where Vq denotes the "cavity Green's function". The 
field ip in Eq. (j35p denotes the Grassmann field '^i in the 
origin, '4> = ^o. 

The cavity Green's function Vq can be determined self- 
consistently from the condition that Eq. (|34p holds for 
any local quantity, including the local dressed propagator 
of the effective lattice theory, 

k 

which, by construction, must coincide with the full 
Green's function of the cavity theory, 

D'°'=V = -{i^i:)s,. (37) 

Then, on one hand, the local propagator on the lattice 
can be determined by using Dyson's equation and thus 
can be expressed in terms of E, as 

D^"" = j2iD"i^r^ -^r^ ■ (38) 

k 

On the other hand, we can also compute the local 
propagator exactly by evaluating the path integral with 
Eq. psp analytically: Using symmetry considerations 
presented in Appendix [Cl we can show that there is a 
U unitary transformation that diagonalizes Vq: 

where d is a 3 x 3 (real) diagonal matrix. We can then 
evaluate the dressed propagator in the local theory by 
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performing a variable transformation ijj ^ (j) = lA^j in 
the intergal, the result being 



D 



loc 



V 



-Do 



l + (.9^-l)\/-detI?o 



(40) 



Furthermore, Dyson's equation also holds for the local 
theory, and moreover in infinite dimensions, the proper 
self-energy E of the lattice theory is the same as in the 
self-energy in the local theory. This follows from the com- 
parison of the skeleton diagrams of the local and lattice 
theories, and the fact that the full Green's functions are 
the same in both theories. We can thus write 



S = I?o 1 - V-^ = -{g^ - 1) v/-detI?o2?(r'- (41) 

With some algebraic manipulations, Eqs. (|38p and ((4T|) 
can be recast in the following self-consistency condition, 




9^ - 1 
V- det S 



J2[D'-\k)-n-K (42) 



This is the central equation of the theory. Having solved 
this equation self-consistently for the self-energy matrix 
S, we can express, e.g., the cavity Green's function as 



-Do 



V-detS ^_i 

0^ — 1 ' 



(43) 



and then determine from that the local Green's func- 
tion using Eq. (j40p or the improper self-energy from the 
Fourier transform of Eq. (P5)l . Remember though that 
the quantities here are 6x6 matrices and inversions mean 
matrix inversions. 

Setting all these together and using the relation 
[Eq. (j30p ] between the density matrix Vij and the im- 
proper self-energy, the kinetic energy can be finally ex- 
pressed as 



while the total double occupancy is 

3 1 



(45) 



a=ll-g2_ /(52_l)^_detE 



This is a good point to discuss the similarities and 
differences with the multiband Gutzwiller method of 
Biincmann et. al^ That approach employs a flexible 
Gutzwiller correlator with multiple variational parame- 
ters. Some of these parameters are allowed to be fixed 
so that the proper self-energy of the effective theory van- 
ishes in infinite dimensions. This leads to important sim- 
plifications when calculating the Gutzwiller expectation 
values. In our formalism, on the other hand, we used a 
single correlation parameter that was motivated by weak 
and strong-coupling theories. Nevertheless, this single 
correlator is able to capture the physics of both limits 
and it therefore accounts for the most important corre- 
lations. Including additional (e.g., two-body correlators) 
should quantitatively modify our results, and change the 
critical value of U , however, we do not expect qualitative 
changes due to them. The advantage of our effective field 
theory approach is that it displays the trionic correlations 
in a rather transparent way, in terms of three-body in- 
teractions. Preliminary results indicate that our cavity 
functional approach could also be modified to use more 
flexible correlators, but this is out of the scope of the 
present paper. 



SOLUTIONS OF THE SELF-CONSISTENCY 
EQUATIONS 

A. Derivation of integral equations 



k a = l 

2 '-,0-1 



1 + .9 
1 



,i?"-i(k)[i?o-i(k)-S]-iS]^ 



{1+9)-- 

where the summation runs only over the first three diago- 
nal elements of the 6x6 matrices. We can also compute 
other local expectation values using the same variable 
transformation [Eq. ([39]) ] as before. The total particle 
density is is given by 

Nt = 3p=i(i^T3V)5. +3(<?2-l)(t)5^ 



a=l i 



-s-^ + 1 



V-dctS 



(44) 



The self-consistency equations are hard to solve in gen- 
eral. However, we can simplify their solution by observ- 
ing that the proper self-energy must have certain symme- 
tries. We therefore need to use only three independent 
parameters to parametrize the self-energy. 



(73 



Si 
El 
Ea 



iO-2 ' 



Eg 

-E3 




(46) 



with (72 and (T3 denoting the second and the third Pauli 
matrices in the Nambu space. Using this matrix, it is 
possible to solve selfconsistency relations (|42)) numeri- 
cally, and we can also calculate the energy. 



E{g, A, P3) A'i2 + K3, + UNDt 



(47) 



where N is the number of lattice sites and we explicitely 
indicated all implicit variables, with respect to which the 
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energy must be optimized. The kinetic energy of the first 
two colors is 

^i2 = 2 2^ek r, , j-2\ ' (48) 



where /k 



(l-Si)(l + /k^) 



is a renormahzed occupation num- 



ber, with 6'k = ^ (ek - Aii2 + \/(ek - Mi2)2 + A^) • The 
kinetic energy of the third color can be expressed as 



[l-S],/(l+g)]% ^, 
= — ^; -f^o , 

i — 2^2 



(49) 



where Kq ~ X]ek</i3 is the noninteracting kinetic en- 
ergy of the atoms in channel 3. Finally, the full double 
occupancy is given by 



2Si + S2 - v/(.92-l)(S^ + E^)S 



l)(S2-t-Si)S2 



(50) 



Notice that the above momentum sums as well as 
the summations included in the self-consistency equa- 
tion contain terms that depend on the momentum only 
through the single particle energy Ck- It is therefore pos- 
sible to turn the momentum sums into energy integrals 
and convert all these expressions into relatively simple 
self-consistent integral equations, involving only the den- 
sity of states 15(e), which becomes the Gaussian on an 
infinite-dimensional cubic lattice, 



D{e) 



1 



(51) 



For high dimensional hypercubic lattices, the qualita- 
tive features of the densities of states are the same as 
in infinite dimensions. In this regard three dimension 
is probably close enough to infinite dimensions, and our 
results should hold away from half filling. One- and two- 
dimensional systems should be investigated by different 
methods since there spatial fluctuations and van Hove 
singularities will play a much more important role. Nev- 
ertheless, even in one dimension, recent Bethe ansata^l 
and density matrix renormalization groupS^ calculations 
seem to support the results presented here and in Ref. ITgI . 



B. Numerical Results 

Our numerical procedure is as follows. First, we fix 
the filling factor p and the interaction strength U . For 
a given set of variational parameters g, A, and p3, we 
use the ansatz (|46)) and solve the selfconsistency rela- 
tions [Eq. (IH])] by iteration32. Once S at hand, we can 
compute the expectation value of the energy i?((7. A, pa) 
using Eqs. (|17|) - ([5(H) . This function can then be mini- 
mized numerically to find the optimum values for g, A, 
and Pi- 

Typical results for a given density pi = P2 = = 1/3 
are shown in Fig. [SI where the energy landscape E{g^ A) 



U=-1.375 t* 



E/|EtI 




-1.02&1. 01 5-1 .005-0.995 
U=-1.625t* 




-1.01 -1.005 -1 -0.995 
U=-1.75 t* 




10 20g 30 40 50 
E/|Et|: ^^^H 



0.9975-0.995 



10 20g 30 40 50 

E/|Ej|: ^^^^^^n 
-0.9995-0.998-0.9965-0.995 



FIG. 5; (Color online) Variational energy E{g,A) compared 
to the energy of the trionic state, Et = SUN p. Only regions 
of E/\Et\ < —0.995 are shown for different values of the 
interaction strength and pi = P2 ~ ps = 1/3. 



is shown as a function of g and A for various values of the 
coupling U. For small values of \U\, the overall energy 
minimum is located at some finite value of the correla- 
tion parameter g > 1 and gap A. This energy minimum 
is shifted to larger and larger values of g as we increase 
\U\. The optimum value of A initially increases with \U\, 
but at larger values of \U\, it starts to decrease again, 
the energy minimum is getting shallower and shallower 
and; eventually, it shifts to g = c» and A = as C/ 
approaches a critical value Uc- This behavior of the op- 
timum values of g and A is similar to the one shown in 
Fig. [51 Within numerical accuracy, the gap vanishes con- 
tinuously at U = Uc, which is characteristic of a second 
order phase transition, i.e., a quantum critical point. The 
correlation parameter g precisely diverges at the same 
value of U, implying that for larger values of \U\ the op- 
timum ground state is a purely trionic state with no su- 
perconducting order. As we discuss below, this scenario 
carries over for all filling factors 0.0333 < p < 0.4833 
that we investigated, where there is always an interac- 
tion strength Uc{p) at which the Gutzwiller parameter 
g diverges. 

In the previous calculation, we fixed all three densities 
to be the same. We have done this in the spirit that in an 
optical lattice, the total number of particles is fixed for 
each color. However, one can decrease the energy of the 
system further in the superconducting phase by letting p^ 
vary and only fixing the total filling fraction p. As shown 
in Fig. [3 the optimal value of the average occupation of 
the third color is slightly less than the filling: p^ < p. 
However, as shown in the inset of Fig. [3 pa approaches 
p as U approaches Uc- This is easy to understand: The 
reason of the formation of an imbalance is that one can 
gain condensation energy by transferring particles from 
color 3 to colors 1 and 2. However, the driving force to 
create this imbalance is the superconducing condensation 
energy. This energy does not coincide with A, but is 
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rather defined as the diS'erence between the the energy 
of the state with A ^ and the energy of the GutzwiUer- 
correlated state with A = 0, 



Econd — E — Eq 



FS 



(52) 



This energy is also shown in Fig. [H] and it vanishes at the 
critical point U = Uc as well. Therefore, the induced 
charge imbalance must also disappear there, in agreement 
with our numerical findings. 

In a system with SU(3) symmetry, this finding implies 
that, although the total number of particles is fixed for 
a given color, one can decrease the energy of the system 
by segregatioufiSi^ i.e., by forming superconducting do- 
mains with unequal numbers of particles and orienting 
the order parameters to point in them in different direc- 
tions, as we already sketched in Fig. [31 Of course, this 
is only possible in a large enough system, where the do- 
main wall energy is compensated by the overall condensa- 
tion energy gain. Interestingly, the tendency to generate 
Pa ^ P can also be viewed as the appearance of a sec- 
ondary ferromagnetic order parameter, 



2" = (cUV) ^ 



(53) 



Therefore we conclude that this superconducting state is 
also a ferromagnet in the SU(3) language. 



C. Breaking the SU(3) symmetry 

So far, we assumed that the scattering lengths are the 
same in all three scattering channels. In most systems, 
however, the three scattering lengths are not equal and 
can vary with an external magnetic field. As we pro- 
posed in Ref. [l^, a possible candidate for realizing the 



G03(p-p3) 
<K> 



0.75 




0.85 0.9 
I U l/lUc I 



FIG. 6: (Color online) Diflterent ground state properties of 
the SU(3) attractive Hubbard model for p = 1/3. For \U\ > 
\Uc\ ~ 1.774i*, the Gutzwiller parameter g diverges, while 
the superfluid order parameter A and the condensation energy 
iScond vanish. Energy is measured in units of t* . 



5: -1.32 




0.39 



FIG. 7: (Color online) One can gain energy by transferring 
particles with color 3 to the superfluid channel: the energy 
minimum occurs for ps < p. The inset shows the optimal 
ps{U) for p — 1. Energy is measured in units of t* . 



trionic state is ®Li. In this system, the magnetic field de- 
pendence of the three scattering lengths has been exper- 
imentally determined.— According to the experimental 
results of Ref. [13, the interaction strengths can be ap- 
proximated in the magnetic field region of 60 - 120 mT 
as 



Up, 



U; 



01 



^01 



B- B, 



01 



(i?-<) 



(54) 

where the parameters i?o, A, and a in this equation can 
be taken from Ref. [13, and U^'^ = Uo, C/q" = 1-23 Uo, 
and = 1.06 Uq. The interaction amplitude Uq = 
Uq^ can be changed by tuning the potential depth. For 
B ^ Br = Bq'^ = 83.41mT, all three interactions become 
approximately the same, and are negative, Uap < 0. On 
the other hand, if we increase the interaction strength 
by approaching the Feshbach resonance at Br, then the 
interaction becomes rather anisotropic in color space. 

Experimentally, there are thus several ways to drive the 
system close to the phase transition regime, i*/|/7o| ~ 1- 
The first possibility is to change the intensity of the laser 
beams, and thereby tune mostly the hopping amplitude 
t* . In this case, one can apply a large magnetic field in 
which the interaction is almost SU(3) symmetrical, al- 
though it is somewhat stronger in channel 12 than in 
the others. The other possibility is to tune the ratio 
t*/|J7o| ~ 1 by changing the magnetic field and approach- 
ing a Feshbach resonance. In this second case, the inter- 
actions can strongly break SU(3) symmetry. 

In both cases, the anisotropic interaction has an impor- 
tant effect: it locks the superconducting order parameter 
so that only A12 ^ 0. The color superfluid phase thus 
becomes a more or less standard U{1) superfluid with an 
additional decoupled Fermi liquid. However, as we show, 
the trionic phase transition can survive in this anisotropic 
limit. 

Let us now investigate the superfluid - trion transition. 
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98 100 102 
B (mT) 

FIG. 8: (Color online) Ground state properties of the Hamil- 
tonian ^ with the interaction strengths (|54p for f*/|L'o| — 
1.43 as a function of B for p = 1/3 fiUing, Uo = Uq'^ . This 
describes a cold atomic system where three hyperfine states 
of ®Li are trapped and the magnetic field is tuned. The inset 
shows the same quantities with the magnetic field fixed at 
B — 116.78mT as a function of f*/|C/o|. Energy is measured 
in units of t* . 



assuming that B > B^- In this region, all interactions are 
attractive, Uap < 0. Perturbation theory tells us that the 
weak-coupling ground state is the color superfluid state^^ 
formed in the channel 12, and in the strong coupling limit 
the ground state is a trionic state. Thus, we can use 
the same type of ansatz as in the SU(3) case. Since the 
structure of the effective theory is uniquely determined 
by the Gutzwiller ansatz state, the effective action, the 
Ward identities and the self-consistency relations for the 
local effective action remain the same as before. The only 
difference is in the expression of the variational energy 
that we need to minimize. 



E{g,A,p3) = K + Y^[Ul2{n^lh^2)G + Ui3{n. 



+ U23{nt2ni3)G] ■ 

The previously used numerical procedure can be modified 
for this case, and we find that the quantum phase transi- 
tion persists even with the anisotropy in the interaction 
strengths. 

The anisotropy in Uap has an important secondary ef- 
fect: it induces somewhat different chemical potentials in 
channels 1 and 2, i.e., an SU(3) "Zeeman field". Deep in 
the superfluid phase, this chemical potential difference is 
not sufficient to create a charge imbalance, pi ^ p2, since 
breaking Cooper pairs requires a finite energyi^ In prin- 
ciple, close to Uc, Ai2 0, and therefore, the chemical 
potential difference could create a first order transition 
to the trionic state. This would invalidate our restric- 
tion pi = p2 in our calculations. However, at the critical 
point trionic correlations also diverge, 5 — > oo, and there- 
fore the effects of this Zeeman field are suppressed. In 
fact, within the Gutzwiller approach, the phase transition 
seems to be continuous even for anisotropic interactions 



(see Fig. [8|). In finite dimensions, on the other hand, we 
cannot exclude a first order transition. 

In Fig. [HI we show different physical quantities com- 
puted for the ground state as a function of B at p = 1/3 
andi*/|L/o| = 1-43. The observed phase transition occurs 
at i? « 94.5mT, and all parameters behave similarly to 
as in the SU(3) symmetrical case. 



VI. ANALOGY WITH QUANTUM 
CHROMODYNAMICS 



The phase transition we found is essentially the 
analogue of the color superconductor - baryon phase 
transition)^ which is believed to occur in QCD. To make 
this analogy with QCD clearer, let us rewrite the parti- 
tion function of the SU(3)-symmetrical Hamiltonian in a 
path integral form, 

S = ^c„(r)^^:^^^„(T)cja(r) 



U 

+ 2 



(55) 



where Q^j^^ — (5,;j(— 9r — Po) — Uj denotes the inverse 
imaginary time Green's function for the non-interacting 
atoms and Uj the hopping between sites i and j. 

We can now decouple the interaction term by using the 
relation 



af3 



directly following from the SU(3) identity, 



(56) 



(57) 



and then performing a Hubbard-Stratonovich transfor- 
mation: 



exp 



= C 




ap 



exp 



at every lattice and imaginary time point. Here we have 
suppressed the imaginary time and color indices of the 
Grassmann variables Ci(r) and the "gluon field", Afi^r), 
and we introduced the coupling 7 = y%5|[/|/2 between 
the gluons and the fermions. 
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Using this transformation wc can thus express the par- 
tition function as 



VA Vc Vc exp 



dr 5<,i[Qa(r),Qa(r),^°(r)] 



In this language, the attractive interaction is mediated 
by real bosonic fields, "gluons", which are ultimately re- 
sponsible for the formation of the trionic ("baryonic") 
phase. Note, however, that in our case, gluons do not 
have a dynamics on their own, rather they are just mas- 
sive Hubbard-Stratonovich fields. As a result, the medi- 
ated interaction is short-ranged, and there is no confine- 
ment. 

In fact, our phase diagram parallels the famous QCD 
phase diagram. In QCD, however, one changes the ratio 
of kinetic energy vs interaction energies through chang- 
ing the chemical potential, while in our case, the hopping 
parameter (i.e., the mass of the Fermions) is changed. 
Therefore, large chemical potentials in the QCD corre- 
spond in our case to large values of t* (i.e., small \U\), 
while small chemical potentials in QCD correspond to 
small t* (large \U\). Another important difference is the 
order of the phase transition: In our case, the two phases 
are separated by a quantum critical point, while in QCD, 
the transition is of first order, due to the long-ranged in- 
teractions. 

Moreover, the superconducting phase wc have is some- 
what different from the one occurring in QCD, since cold 
atoms in optical lattices have only a color quantum num- 
ber, while quarks in QCD have additional flavor degrees 
of freedom. The color superfluid emerging in our case has 
a nontrivial SU(3) color content and is analogous to the 
color superconducting phase in two-fiavor QCD where 
only two flavors of light quarks are considered ^ In 
the alternative theoretical scenario of three-flavor QCD, 
the superconducting state is expected to be color-flavor 
lockedi^ 



VII. CONCLUSIONS 

In this paper, we investigated a system of fermionic 
cold atoms with three internal quantum states and at- 
tractive interactions among them. We have shown by 
heuristic arguments as well as by a detailed variational 
calculation in the limit of d — > oo-dimensions that this 
system shows a quantum phase transition from a cor- 
related superconducting state to a trionic state, where 
fermions form three-body bound states and the super- 
conducting order disappears. 

We also showed that in the case of SU(3) symmetry, 
ferromagnetic ordering appears as a secondary order pa- 
rameter. In a closed system of atoms, where the number 
of atoms is conserved for every color, this implies seg- 
regation, i.e., the formation of domains where the color 



of the atoms forming the superconducting state changes 
from domain to domain and the density of the atoms 
involved in the superconducting order is slightly higher 
than that of the unpaired ones. 

We have also demonstrated that the phase transition is 
robust against breaking the SU(3) symmetry of the inter- 
action. This symmetry breaking influences the structure 
of the superconducting state in that it pins the order pa- 
rameter to the most strongly interacting channels, but 
docs not influence the existence of the phase transition 
itself. 

We have not studied, however, the case where an un- 
equal number of atoms is loaded for the three hyperfine 
states. This situation has been experimentally studied 
for a two-component ^Li system.— There, a segregated 
superconducting state has been observed in the center 
of the trap, where the two densities become equal in or- 
der to form Cooper pairs, while at the edge of the trap, 
spin densities are different. A similar scenario is expected 
for a slight imbalance of the three densities in our case. 
It is, however, also possible that for SU(3) symmetrical 
interactions, three different superconducting domains of 
unequal size can form for large enough traps. Forma- 
tion of a Fulde-Ferrel-Larkin-Ovchinnikov state cannot 
be excluded either, although this state has not been ex- 
perimentally observed 

Here, we should also remark that we solved the 
Gutzwiller problem only in d = oo dimensions, where 
it predicts a phase transition. Although the phase tran- 
sition survives in finite dimensions too, there we do not 
know if the variational Gutzwiller wave function displays 
a phase transition or not. As long as g remains finite, one 
has a step in the momentum distribution at the Fermi 
energy and the system can gain energy by forming a 
superconductor. Therefore, to have a phase transition, 
the parameter g must diverge at a finite critical value 
of U . This certainly happens within the local approxi- 
mation (Gutzwiller approximation) in finite dimensions 
too, and one can argue that this should also happen if 
one could evaluate the finite-dimensional Gutzwiller ex- 
pectation values exactly. Unfortunately, restricted Monte 
Carlo calculations would be needed to give a definite an- 
swer to this question. 

Finally, let us briefiy discuss how the transition from 
a superfluid to a trionic state could be experimentally 
detected. One way to observe the condensate is by de- 
tecting vortices. In a perfectly SU(3)-symmetrical system 
there are no vortices, because any vortex can be twisted 
away due to the large internal symmetry of the order 
parameter However, if one breaks the SU(3) symmetry 
down to U(l) then the superconducting phase will con- 
tain usual vortices that can be relatively easily observed 
optically;^ In case of ®Li, e.g., this can be achieved by 
approaching the Feshbach resonance at B^- = 83.41 niT 
from the high-field side, and thereby having a condensate 
in channel 12. Then our prediction is that the superfiuid 
state and thus the vortices disappear for large enough 
\Uafj\l't*- This ratio can be tuned by cither changing 
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the laser beam intensities and thereby the amplitude of 
the optical lattice or by changing the external magnetic 
field. A further consequence of having U13 ^ U23 is that 
the interaction induces a chemical potential difference for 
channels 1 and 2. Although we do not see a sign of it 
within the infinite-dimensional Gutzwiller approach, this 
may well render the phase transition first order in finite 
dimensions. 

Pinning the order parameter to a single channel has a 
further advantage. Since one knows that the order pa- 
rameter is simply in channel A12, one can use the method 
of Ref. [3^ to measure A12 by sweeping the magnetic field 
through the Feshbach resonance at Br = 83.41mT and 
detecting the momentum distribution of the molecules 
thus formed. This method would enable one to detect 
|A| as a function of lUapl/t* and verify the existence of 
a critical value where it goes to zero. 

Stability of the three-color ^Li system may also be an 
important issue. High magnetic fields will probably con- 
siderably stabilize the three-color condensate, but using 
mixtures of different fermionic atoms with all scattering 
lengths being negative may also be a possibility. In such 
a composite system, other interesting phenomena could 
also take place due to the difference in atomic masses. 
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APPENDIX A: WARD IDENTITIES 

Here, we connect certain expectation values with the 
self energy S defined by Eq. (|24p to express the kinetic en- 
ergy in terms of the self-energy itself. First, we shift the 
integration variables of F in Eq. ((22|) following Eq. (|26l). 
This transformation automatically transforms the conju- 



gate fields by 

^ = + A^/pTaT^o,. (Al) 
p 

Using these, the linear and quadratic terms become 

= ^ E ^^^i^'^j + (1 - ^) E ^^^3$, 

+ ^^y^E^»^3^?;^3/, ■ (A2) 

Transforming the interaction term is more complicated. 
First, we transform the on-site density, 

*»r3*^ = $,r3$, -2Al>,r3^i?°„T3/^ 

m 

+ A2^/„T3Z?°,r3i?LT3/m , (A3) 

nm 

where we used that ^iT^D^jT^Ij = IjT^D^^^Tz^i. Here, 
we note that all terms behave like scalars with respect 
to both the Grassmann algebra and the color matrices; 
and they commute with each other. When performing 
the functional derivation, we will have to expand the ex- 
ponential in terms of / to the order 0{P). Thus, when 
calculating the third power of Eq. (|A3[) . we can drop cer- 
tain higher order terms. The interaction term can then 
be expanded as follows: 

E i^rT^^rf 

r 

« ^($,T3$.)3 - 6A^($,T3$,.)2($,.r3A/->3/m) 

r r.ra 

+ 12A2 (/„r3M-;r3$r)(^rr3$.)($.r3Af-ir3/m) 

r,n,m 

+ 3A2 («'rr3$.)'(/„r3M„-V3il/-^r3/„) + • • • . 

r,n,m 

Now, we can expand the exponential in the definition 
of the generating functional. Introducing the coupling 
u = — 1 oi the effective theory, we obtain the following 
long expression 



F = In J Vf]V7] e-^+^^''^''^' « 

// A^ — 2A - 1 - 

VaVa e--^ I 1 + E -^'■^sAC' ^s^r' + ^\^u ^ {Ir'TzM^.lr^n^ $,)($,T3M-iT3/." ) 
\ rr' j-^'^" 

+ \\^u Y {Ir'r3M-,ld^T3M~;,r,I,.^) + hi - \f ^(/,r3$,,)(|.,,r3/,0 



13 



rr'pp' 

-iA(l - \)U iIrT3^r)id^,^r'T3M-,l.T3lr'^) - ^Xil - X)u ^ (/.rgM" } Tgd,^ <I>.0 " Tg/." ) J , 

rr'r" ^-p'^" / 



(A4) 



where a, a are the transformed Grassmann fields, which correspond to $. Now one can perform the fmictional 
derivation to get the dressed propagator, 



SJiSJj SJiSJj 
SJiSJj SJ.SJj 



T3 



5151 



T3 = (2A - \^)Dl - (1 - A)2($,#^)5 - \^uY,DlM<^r^r)sT3Dl^ 



r rr' 

+ (1 - A)Aw^($,<$,)5r3Z?,"^ + (1 - \)\uYD°M<^rd^<tj)s • 



(A5) 



Using that D = — we can rewrite this in the foUowing form: 

A 



A 



A, = Dl-^-^uJ2DiMn^<^r^r)sr3D°r,-l^uYDtr3{d^)sD[ 
A 



2- A 



-u 



2 -A" 

r 

Y i^°.T3($,df d^,^r')sT3D°,^ 



+ ^uY{'^>.d^^r)sr3D"r, 



L-^uYDtT,{<i>,d^^,)s. 



(A6) 



Substituting A = 1 leads to 



-M^I3°,r3(nf$,$,)5T3Z?°, + 

r 

-U^Z?^,T3(df)5A",+ 

r 

-1/2 Y DtT3{<^>rd^d^,^,.')sr3D"r'j, (A7) 



while setting A = we obtain another interesting iden- 
tity, 



+iE^-^3(<i>.<*,)5. (A8) 

r 

Compared these to the definition of the self-energy 
Eq. jSll), one finds the identities Eqs. ^ - ^ and 
thus the expression of the kinetic energy in terms of the 
self-energy. 



APPENDIX B: DERIVATION OF THE 
FUNCTIONAL FORM OF THE LOCAL ACTION 



The local action is defined by 



1 



-Sl 



(Bl) 



where the prime means integration over all Grassmann 
variables except for those defined at the origin. We shall 
derive its functional form in = oo dimensions. It is 
useful to collect the terms in the action into three groups, 

S = Sq + S' + S^°\ 



where 



5o 



(B2) 



*o[^""']oo*o-(,9'-l)io (B3) 



5(0) 



2 ^ 



^AD'~%-^,^{9'--l)yU (B4) 



contains only terms at the origin, 

1=^0 

is the effective action on a lattice without the origin, and 
S' = -ij] [*4^""'].o*o + *o[i5"-']o,*,] (B5) 
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are the terms which connect the lattice and the cavity. S' 
can be viewed as a generating functional term with the 
currents hi = ^o[D^~^]oi and their adjungates. Thus, 



In J V'fjV'ij e 



n ii...i„,ji...j„^a 



(B6) 



since only connected graphs arc generated this way. Let 
us first restrict ourselves to sites iq and jq being only 
nearest neighbors, and consider terms of order n. If all 
site indices are different, then summation gives a factor 
d^". Since each bare propagator D""^ is proportional to 
t ~ l/Vd, they give a prefactor ~ c?~". The correlation 
function is connected, and thus the distance between ex- 
ternal sites is I |i — il I > 2, and the scaling of the hopping 
ensures that the correlation function is 



n 



\ conn 



< const/d2"-i. (B7) 
terms give a contribution of the order 



Thus, these 
l/d"-\ 

If two of the labels {iq} and {jq} are the same, then 
the summation gives a factor d^""^, however, connected 
diagrams scale at most l/d^"~^, and the final result is 
again ^ l/d"^^. This argument can be generalized to 
any combination of identical labels. This means that 
only the quadratic terms n = 1 survive the d — *■ oo limit, 
which can be added to the quadratic term in Sq, leading 
to the bare propagator of the local theory Vq- The self- 
consistency relations are required to determine its value. 

This line of argumentation can also be extended to the 
case where iq and jq run over non-nearest neighbor sites. 
One only has to use the property that (1?°)^^ falls off 

faster than l/c?ll*~^ll/^, with ||i — j|| the number of steps 
needed to reach the lattice site j from site i. 



APPENDIX C: DIAGONALIZATION AND 
SYMMETRY PROPERTIES OF THE LOCAL 
ACTION 



This structure of V^^ enables us to say something about 
its eigenvalues. Let us assume that Vj is a right-hand 
side eigenvector of V^^ (a column vector). 



(C3) 



Then, Wj = tiV* is also an eigenvector, 

(C4) 

since the eigenvalues are real. As a consequence, we can 
construct a unitary matrix U'^ as 



(vi , V2, V3, Tiv*, riv^, riv*) 



which transforms Dq to a diagonal form, 



(C5) 



(C6) 



with d is a real diagonal matrix. 

Now, we can show or see that U is a well-defined trans- 
formation for ^ in the sense that the identity ^ = ^^ti 
also holds for the transformed eigenvectors, $ ^U'^,^ = 
^U+, and 



(C7) 



Note that this identity is necessary to carry out the path 
intergral in terms of the components of the transformed 
field $ as new independent variables. This equality is 
satisfied if 

= $ = $^Ti = (ZY«')'^Ti = ^-^Z^^ri, (C8) 
i.e., if the matrix U satisfies the condition 



(C9) 



Equation (jC9[) can readily be verified by using the prop- 
erties of the eigenvectors 



In this Appendix, we show that the matrix U that 
diagonalizes the local propagator "Dq^ has a special sym- 
plectic symmetry, which implies that one can evaluate 
the Green's functions by performing a variable trans- 
formation with this matrix. We first observe that the 
Nambu spinor satisfies the relation ^ = '^^ti where 
Ti = ui ® Sap and tri is the first Pauli matrix acting in 
the Nambu space. Therefore, we can parametrize in 
the following way. 



^-1 _ 



H 

A+ 



A 



(CI) 



where H is Hermitian and A is antisymmetric. Thus, the 
following identity holds 



(C2) 



TlU'^Tl ^ Tl (V J , , V5 , Tl , Tl V2 , Tl V3 )ri = 

= n (n , Tl V2 , Tl V3 , V* , , ) 

= (Xl , V2 , V3 , Tl V J , Tl , Tl V3 ) = 

= li+ . (CIO) 

Relations [Eqs. (|C6P and (jC7P ] together with these imply 
that the Green's function of the local Green's function 
can be computed by first transforming to the diagonal 
basis of and then transforming back the Green's 
function to the original basis by using the transformation 
U. Importantly, the interaction term of the local action is 
invariant under the transformation \1/ — > This follows 
simply from the facts that the interaction can also be 
expressed as ^ 11^1=1 ^^'^ ^^^^ determinant of U 
is just 1. 
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